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ABSTRACT 

A simple model of the coronal magnetic field prior to the CME eruption on May 12 1997 is de- 
veloped. First, the magnetic field is constructed by superimposing a large-scale background field 
and a localized bipolar field to model the active region (AR) in the current-free approximation. 
The background field is determined from the normal component of the observed photospheric 
magnetic field averaged over the longitude of the Sun. The AR field is modeled with the help of a 
subphotospheric dipole whose strength, location, and orientation are optimized to fit the magnetic 
field obtained from an MDI magnetogram. Second, this potential configuration is quasi-statically 
sheared by photospheric vortex motions applied to two flux concentrations of the AR. Third, the 
resulting force-free field is then evolved by canceling the photospheric magnetic flux with the 
help of an appropriate tangential electric field applied to the central part of the AR. 

To understand the structure of the modeled configuration, we use the field line mapping 
technique by generalizing it to spherical geometry. It is demonstrated that the initial potential 
configuration contains a hyperbolic flux tube (HFT) which is a union of two intersecting quasi- 
separatrix layers. This HFT provides a partition of the closed magnetic flux between the AR and 
the global solar magnetic field. Such a partition is approximate since the entire flux distribution is 
perfectly continuous. The vortex motions applied to the AR interlock the field lines in the coronal 
volume to form additionally two new HFTs pinched into thin current layers. Reconnection in 
these current layers helps to redistribute the magnetic flux and current within the AR in the 
flux-cancellation phase. In this phase, a magnetic flux rope is formed together with a bald patch 
separatrix surface wrapping around the rope. Other important implications of the identified 
structural features of the modeled configuration are also discussed. 

Subject headings: Sun: coronal mass ejections (CMEs) — Sun: flares — Sun: magnetic flelds 



1. INTRODUCTION 

The May 12 1997 CME occurred during solar 
minimum in an isolated active region near disk 
center. Since the solar magnetic field during this 
period had minimum complexity, this event is es- 
pecially favorable for detailed study. Therefore, it 
was selected by several working groups (SHINE, 
CISM and MURI) for modeling and analysis. Sev- 
eral observational features of the event have been 
reported (Hudson et al. 1998; Plunkett et al. 1998; 
Sterling et al. 2001; Thompson et al. 1998; Webb 



et al. 2000; Gopalswamy & Kaiser 2002; Lundquist 
et al. 2003; Arge et al. 2004; Li et al. 2004; Liu 
2004; AttriU et al. 2006, 2007), which has provided 
a good base for creating a magnetohydrodynamic 
(MHD) model of the event. For the latter, we have 
used the following three-step approach. 

First, a potential magnetic field is determined 
from the observed photospheric distribution of the 
normal field component. Second, the configura- 
tion is sheared quasi-statically by photospheric 
vortex motions applied to the opposite polarities 
of the active region (AR). Third, an appropriate 



tangential electric field is applied then to the cen- 
tral part of the AR to cancel half of the photo- 
spheric magnetic flux. This third step is moti- 
vated by the observed flux-canceling properties of 
the photospheric motions in the event (Li et al. 
2004; Ambastha & Mathew 2001). Theoretically, 
the third step is intended first to form a fiux rope 
and then to destabilize it by detaching the rope 
field from the photospheric one, which has even- 
tually to produce its eruption into the corona and 
heliosphere. This complete modeling program has 
begun by our group and it will be described in 
forthcoming papers. Preliminary results for these 
simulations have been presented by Linker et al. 
(2006). In this first paper, we consider a simpli- 
fied model of the magnetic field and its structure 
at each of the above three steps in order to iden- 
tify the generic structural features of the configu- 
ration prior to its eruption. By the term "generic" 
we mean those features of the structure that are 
likely to remain in more accurate models of the 
AR magnetic field (AR8038 at N21W08). 

Observationally, this field had an apparently 
simple bipolar structure (Liu 2004) dominating 
over a more complicated but distant and weak 
surrounding field. As a starting point, it is nu- 
merically extrapolated in the potential approxi- 
mation from the synthetic magnetogram obtained 
by merging an MDI magnetogram with the cor- 
responding synoptic map. Prom the methodolog- 
ical point of view, it is sensible, however, to be- 
gin with a simpler model by approximating the 
numerical potential field by a sum of two compo- 
nents, an azimuthally averaged large-scale mag- 
netic field plus a bipolar field of the active re- 
gion itself. This model, comprising essentially two 
bipolar field components with different scales and 
orientations, has a deceptive simplicity. It turns 
out that the structure of such a field is nontrivial 
even at the initial current-free state, and it be- 
comes even more complicated at the shearing and 
flux-canceling phases, features that have impor- 
tant implications for understanding the event. 

The analysis of this structure prior to eruption 
is a major goal of the present paper. The paper 
is organized in the following way: the construc- 
tion of the field model is described in sections 2 
and 3; in section 4, the so-called squashing factor 
is generalized to the case of spherical geometry 
in order to investigate in section 5 the structure 



of the modeled field. Section 6 describes the re- 
lationship between field line connectivity and the 
current distribution. Our results are summarized 
in section 7. Some details of the initial potential 
field model are provided in the Appendix. 

2. THE INITIAL POTENTIAL MAG- 
NETIC FIELD 

Although the computation of potential mag- 
netic fields is routine, the investigation of the field 
structure and its impact on the subsequent MHD 
evolution is a novel approach. For this reason, we 
prefer to idealize our initial field Bq as much as 
possible by keeping only the most essential prop- 
erties of the potential field computed from mag- 
netogram data. Since our primary interest lies in 
exploring the eruption process, the initial field Bq 
ought to incorporate at the very least the large- 
scale coronal field Bq, which we call the back- 
ground field, and the field Bar of the active re- 
gion itself. Thus, the resulting field can be written 
as 

Bo = Bq + Bar = -VFq - VFar, (1) 

where Fq and Far are harmonic potentials of the 
magnetic field. It should be noted that although 
both fields have a similar bipolar character, then- 
have in general different orientations and length 
scales: the background field varies on the scale of 
the solar radius Rq, which greatly exceeds the size 
of the active region L. 

2.1. Background field Bq 

For the sake of simplicity, we also assiime that 
^0 and, in particular, its photospheric radial 
component -Br0lr=fl„i is axisymmetric. The lat- 
ter can be approximated by averaging the ob- 
served Br{RQ,0, 0) over the longitude 4> and then 
by smoothing and fitting it to the sum 

JV 

B.0U^^= 5^ c„cos"^, (2) 

n— no 

where the coefficients c„ are determined by a least 
square minimization. Imposing additionally the 
flux-conservation condition and sacrificing the ac- 
curacy of fitting at the poles, where the observa- 
tional data are not reliable because of the projec- 
tion effect, we have found that equation (2) pro- 
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vides at no = 5 and = 11 quite a good approx- 
imation to the large-scale field of the Sun (Fig- 
ure 1). 
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Fig. 1. — The latitudinal distribution of the photo- 
spheric radial magnetic field after averaging over 
longitude, together with its smoothed analytical 
approximation (thin solid and thick grey-shaded 
lines, respectively). 

Under these assumptions, the potential of the 
background field can be expressed in terms of 
spherical functions as 



+ C+r")P„(cos^), (3) 



where P„ (cos 0) is a Legendre polynomial of n-th 
order. The terms with n = are omitted here 
because the total magnetic flux must be zero and 
Fq is deflned up to an arbitrary constant. 

The unknown coefficients in this expression 
are determined by linear algebraic equations de- 
rived from the boundary or asymptotic conditions. 
The first number iV of such equations can be ob- 
tained by equating the corresponding linear com- 
binations of C^, k = 1, . . . , n, at cos" 6 in 



B, 



dFr. 



dr 



(4) 



r=Rc 



to either or c„ depending on n (see equation 
(2)): the first and second case corresponds to 



n — 1, . . . , TT-o — 1 and n — hq, . . . , N, respectively. 
These equations are rather cumbersome, but they 
can easily be retrieved, so we do not present them 
here explicitly. 

The form of the remaining iV equations for 
depends on what kind of behavior at large radii 
r we would like to have for JBq. The simplest 
choice is to require Bq\^_^^ — > 0, for which the 
corresponding coefficient equations are simply 



C+=0, n = l,...,N, 



(5) 



- and this is one of the two cases that we have 
considered. 

However, for modeling the solar wind effect, the 
second choice of boundary condition, namely. 







with the source-surface radius Rgs equal, say, to 
2.5 -Rq, is also of interest. According to equation 
(3), this yields 

R-^-'C-+RlC+ = 0, n = l,...,N, (6) 

which reduces to equation (5) at i?ss oo, as ex- 
pected. The magnetic field under this condition 
becomes radial at the sphere of radius Rss, and 
the corresponding approach is called the potential 
field source surface (PFSS) model (Altschuler & 
Newkirk 1969; Schatten et al. 1969). This type of 
solution models an open field structure sustained 
by the solar wind, which makes it possible to esti- 
mate its impact on the coupling of Bq and Bar- 
For the analysis of the field line structure in 
both configurations, it is useful also to have an 
explicit expression of the magnetic flux function 
^'©(r, 9), which in our case is 



N 



*0 = E 



n=l 







[(- 






n J 



[cos^P„(cos6i)-P„+i(cos6i)]. (7) 

The contours of \t'0 in the planes </> = const, define 
the field lines of Bq, since 



Br 



V*0 X 



(8) 



We chose Rq to be a unit length scale in the 
problem, so that the coefficients of both type of 
solutions are calculated from the above equations 
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at R-., = 1. Table 1 in appendix A presents their 
numerical values and shows that the coefScients 
C~ of these two solutions are not very different, 
and the coefficients of the second solution arc 
rather small. This suggests that the solar wind 
does not much affect the magnetic field in the 
lower part of the corona, in agreement with our 
study of the entire magnetic field Bo described 
below. It is also in agreement with the structure 
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Fig. 2. — Magnetic field lines of the axisymmetric 
background potential field Bq in the plane a; = 
at r < Rss — 2.5i?0 for two different boundary 
conditions (a) BqI^^ and (b) -FqI^^h^^ = 
0. The null points of the second configuration are 
indicated by Y letters turned by ±90°. 

of the field lines plotted as equidistant contours 
of the flux function in the plane a; = of a 
Cartesian system of coordinates {x,y,z). Figure 



2 clearly demonstrates that both structures are 
nearly the same everywhere accept for r > l.Si?© 
at equatorial latitudes. As expected, the source- 
surface configuration has a null line along the 
equator of the sphere r = Rgg to accommodate 
the direction reversal of the open field. 

2.2. Active-region field Bar 

The normal component of the photospheric 
magnetic field in our AR has a well-pronounced 
bipolar character. It is natural therefore to ap- 
proximate the corresponding potential field with 
the help of a fictitious subphotospheric dipole, 
whose potential is 



This analytical expression was used to fit the nu- 
merical potential field solution B^um computed 
from the original magnetogram data. We have 
optimized the location of the dipole and its 
moment m-d to minimize the difference between 
BrQ + Br AR and Br num. More precisely, ra and 
rrid are fitted to minimize the sum of (B,-© -|- 
Br AR— Br num)^ evaluated at 14 points which were 
randomly scattered over the AR on the sphere 
r = Topt- The radius Topt is arbitrary to a certain 
extent, and its value was chosen equal to 1.12i?0 
to filter out magnetic inhomogeneities of the size 
< O.I-R0 in the lower part of the solar atmosphere. 

In spite of its simplicity, the described analyt- 
ical model provides a relatively good approxima- 
tion of the numerical potential field at r = Vopi 
(see Figure 3). The model becomes even more ac- 
curate at larger r, because higher-order harmonics 
d(^sc,ribing the small-scale inhomogeneities in the 
numerical solution decay faster with r than the 
lower-order ones. 

Table 2 in appendix A presents the resulting op- 
timized parameters of the dipole in spherical co- 
ordinates for both types of the above boundary 
conditions. It shows that the magnetic moments 
rrir of these two types of solutions differ less than 
2.4% from each other, while the differences be- 
tween other parameters are even smaller. This is 
consistent with the above conclusion that the so- 
lar wind does not greatly influence the magnetic 
structure in the low corona. 
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3. QUASI-STATIC MHD EVOLUTION 
OF THE MAGNETIC CONFIGURA- 
TION 
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Fig. 3. — Contours of Br ar and Br num. — BrQ 
(grey-shaded) at r = ropt in the ^)-plane for 
two different boundary conditions (a) Bq\^^^ — > 
Oand (b) FqI^^^^=Q. 



The above potential field with the boundary 
condition Bq\^^^ is used as the initial con- 
figuration for our simplified MHD model of the AR 
evolution prior to eruption. This field is not in- 
tended to describe the coronal field at any instant 
of its evolution, but rather to serve as a convenient 
starting point for further modeling. To energize 
this configuration, it is quasi-statically sheared by 
the photospheric flows that preserve the initial 
photospheric radial component of magnetic field. 
This step produces a stressed configuration with 
a nearly force-free magnetic field. Its subsequent 
evolution is driven by imposing a special photo- 
spheric distribution of the tangential electric field 
that models the observed flux-cancellation pro- 
cess. 



3.1. Shearing/twisting 
flux spots 



of the magnetic 



The initial potential field Bq is sheared by pho- 
tospheric flows such that its footpoints are rotated 
in the flux spots of the AR in opposite directions 
by following approximately the contours of the ini- 
tial photospheric distribution of the radial field 
component Borp,{4',d) = Borlr^R^- "^^^ corre- 
sponding shearing flow is described by 
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B: 



-Vt X 



Ore 



(10) 



where Vt stands for the gradient operator tangen- 
tial to the photosphere. This velocity field van- 
ishes at the PIL defined by BorQ{<l>,0) = and 
preserves the initial distribution of Bor^, since 

Vf (Bor^Vsh) = 0. 



The function 



/ = a{t) exp 



sm 7 cos 7 



+ 



(11) 



provides an oval mask concentrating the shear in 
the AR near the PIL, where a{t) is an ampli- 
tude linearly ramping until a{t) = 0.025 Rq/t is 
reached at t = 0.5 r. After that a{t) remains con- 
stant until t = 3.0 T, when it drops linearly to zero 
att = 3.05 r and remains zero from there on. The 
time unit r must be large enough to guarantee 
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quasi-static evolution of the whole magnetic con- 
figuration. Specifically, it was chosen to be ~ 1800 
times the Alfven time in the AR, and 4 times the 
Alfven time in the polar field region. 

The angle variables are measured hereafter in 
radians and defined as 



7 
5 



7o -|- arctan (f?o 
R-Rq 



= arccos ■ 



Rq 



>o), (12) 
(13) 



where 70 ~ —1.396 is a tilt of the oval, R and Rq 

are the radius-vectors of the photosphcric points 
with longitudes and 0o = 3.026 and latitudes 
6 and ^0 = 1-184, respectively. The parameters 
W0 and controlling the sizes of the oval mask 
are chosen to be we = 0.06 and = 0.02. The 
resulting photospheric velocity field determined by 
equations (10) (13) is shown for the described set 
of parameters in Figure 4a. 

3.2. Model of the flux cancellation 

The process of flux cancellation is modeled by 
imposing at the photosphere the following tangen- 
tial electric field 



(14) 



where the flux function if} is determined from the 
two-dimensional Poisson equation on the sphere 
r = Rq such that 



Or©- 



(15) 



The choice of a spatially uniform factor g provides 
a simple uniform reduction of -Bor© - However, the 
flux function ij) and the related electric field 
decay in this case too slowly with distance from 
the AR compared to the respective magnetic field. 
This leads unfortunately to the undesirable forma- 
tion of a strong boundary layer covering a substan- 
tial area outside of the AR. To localize within 
the AR, the mask function 



g = c{t) exp 



(16) 



is used with = 3.022 and ^* = 1.187 repre- 
senting an approximate position of the AR cen- 
ter, and w = 0.065 determining the extent of the 
electric-field localization. The amplitude c{t) here 
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Fig. 4. — (a) The vector field of the photo- 
spheric velocity in the shear/twist phase and the 
respective contours of i?r|r=i?Q = 44 i G, i = 
—5, . - - , 6; solid and dashed lines correspond to 
positive and negative values, respectively, while 
the thick dashed line represents the PIL. (b) The 
longitudinal profiles of Br\^_j^^ in the AR at the 
middle latitude 6^, = 1.187 for five equidistant 
time moments: the darkness of the grey-colored 
curves is monotonously increased from the initial 
to the final profile. 
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is linearly increased from t = 3.1 t until some 
c{t) = Cmax is reached a.tt = 3.3 t. After that, c{t) 
remains constant until t = 7.5 r, when it drops lin- 
early to zero at t 7.7 r and remains zero from 
there on. The constant Cmax is chosen from the 
requirement to reduce the AR flux approximately 
to a half of its initial value in the specified inter- 
val of time. The resulting reduction of the pho- 
tospheric normal magnetic field in the described 
flux-cancellation model is depicted in Figure 4b. 

4. FIELD LINE CONNECTIVITY IN 
SPHERICAL GEOMETRY 

As we will see, our initial configuration has a 
simple topology such that no magnetic null points 
(NPs) or bald patches (BPs) are present in the 
modeling AR. BPs are segments of the PIL at 
which field lines are concave towards the corona 
(Titov et al. 1993). An exception is the above- 
mentioned null line at the source surface r = i?ss, 
which accommodates the direction reversal of the 
open field at the equator in the PFSS model. 
This null line, however, is outside of the AR, so 
its structure may have only geometrical features 
relevant for understanding the eruption process 
under study. The features of major interest are 
quasi-separatrix layers (QSLs, Priest & Demoulin 
(1995)) and, particularly their combination in the 
form of two self-intersecting QSLs, the so-called 
hyperbolic flux tubes (HFTs, Titov et al. (2002)). 
This is because HFTs have properties favorable for 
developing current sheets and reconnection (Titov 
et al. 2003; Galsgaard et al. 2003; Aulanier et al. 
2005), even if the respective configuration is free 
of NPs and BPs. The properties of such struc- 
tural features and their manifestation in observed 
events are recently reviewed in detail by Demoulin 
(2005, 2006). 

These features are determined with the help of a 
point-wise analysis of the magnetic field line map- 
ping. For small ARs, where the photosphere can 
be approximat(^d by a plane, it is convenient to de- 
fine this mapping in a Cartesian system of coordi- 
nates (x, y). Such coordinates allow us not only to 
determine the location of the footpoints but also to 
measure the distances between them in the x and 
y directions. The mapping {X{x,y),Y{x,y)) de- 
fin(^s the connection between one footpoint {x, y) 
and another {X, Y) by a magnetic field line. No- 



tice that this mapping is defined only for closed 
field lines, while its generalization for open field 
lines requires an introduction of an upper bound- 
ary surface (sec below). The local properties of the 
field line mapping are described by the Jacobian 
matrix 

(17, 

whose elements, however, are not invariant with 
respect to the direction of mapping. Therefore, 
they cannot be used directly for characterizing 
the field line connectivity and for detecting QSLs 
in configurations. One can prove that a natural 
quantity which satisfies this requirement in the 
case of a plane photospheric boundary is expressed 
in terms of the elements of D as (Titov et al. 2002) 

Q = A^V|A|, (18) 

where 

Ar2 = a2 + 62^c2+d2 (19) 

is the so-called norm squared (Priest & Demoulin 
1995) and 

A = ad -he (20) 

is the Jacobian of the mapping. To explain the ori- 
gin of Q, note that the Jacobian matrix represents 
a linearized field line mapping in the vicinity of a 
given footpoint. Therefore, any infinitesimal cir- 
cle centered at such a point will be mapped along 
the neighboring field lines into an infinitesimal el- 
lipse at the other footpoint. It can be shown that 
the value Q/2 + \JQ'^ — 1 determines the aspect 
ratio of such an ellipse [see the details in (Titov 
2007)]. It characterizes the degree of squashing 
of the corresponding infinitesimal magnetic flux 
tubes and turns simply into Q for Q » 2. Since 
QSLs are determined by large values of Q, we will 
neglect the slight difference between Q and the 
above exact value of squashing degree by calling 
Q itself the squashing degree or squashing factor 
(Titov et al. 2002). 

Another independent geometrical quantity in- 
variant to the direction of mapping (up to the sign) 
is even simpler than Q. It determines the degree of 
expansion or contraction of infinitesimal flux tubes 
and, due to the flux conservation in them, can be 
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expressed in terms of the ratio of the normal field 
components at the connected footpoints as (Titov 
et al. 2002) 



or at the source surface, respectively. Thus, for 
the elements of matrix (17), we obtain in our case 



i^ = iog|A| = iog|s„/i3;; 



(21) 



Here and Bn are normal components of the 
magnetic field at the conjugate "target" and 
"launch" footpoints, {X{x,y),Y{x,y)) and {x,y), 
respectively. In practice, the numerical calcula- 
tion of A through this ratio is more precise than 
that given by equation (20) and therefore it should 
be used for computing Q in equation (18) as well. 

The quantities Q and K provide a valuable tool 
for investigating the geometrical structure of coro- 
nal magnetic fields. For example, they made it 
possible to convincingly demonstrate that HFTs 
arc a generalization of separator field lines to the 
case of topologically simple magnetic configura- 
tions having no NPs or BPs (Titov et al. 2002). 
This is particularly interesting for our study, since 
HFTs appear naturally in quadrupole configura- 
tions (Titov et al. 2002), which can be thought 
of as being the superposition of two nested dipolc 
fields — exactly the type of potential field consid- 
ered in the previous section. 

Such a field-line mapping technique can be gen- 
eralized in covariant form to the case where the 
magnetic configuration is bounded by curved sur- 
faces of arbitrary shape (Titov 2007). Here we 
will restrict ourselves to spherical boundaries and 
use a much simpler heuristic approach to derive 
Q in a global spherical system of coordinates. 
Notice first that our coronal magnetic field fills 
the volume between two spheres of radii Rq and 
Rss- For measuring distances between infinites- 
imally close footpoints, one needs to introduce 
two-dimensional Euclidean bases, which are tan- 
gential to these boundaries at each of their points. 
In the traditional spherical system of coordinates 
with the latitude ^ = or tt at the poles of 
the Sun, these must be pairs of orthoiiormal vec- 
tors {e^,e$). Then the infinitesimal increments 
of length along and eg in one of the two po- 
larities are 5x = Rq sin 6 d(p and Sy = Rq d9, 
respectively. Let {^{(j),6),Q{<j),9)) be the vector 
function representing the field line mapping, then 
similar increments at the conjugate footpoint are 
5X = i?* sin 9 d<I> and SY = i?, d6 with either 
i?* = Rq or i?* = Rss depending on whether the 
conjugate footpoint is located at the photosphere 



a = 



R^ sin e 9$ 
Rq sin 9 d(j) ' 



c = 



R* dQ 
Rq sin 9 d(f) ' 



d = 



Rq 89 
R^de 
Rq 39 ■ 



(23) 



It is these expressions that must be used for com- 
puting norm (19) in squashing factor (18) general- 
ized for the case of both open and closed magnetic 
fields in the coronal volume restricted by spheri- 
cal boundaries. The value of A as a part of the 
whole expression for Q is more accurately com- 
puted through the ratio of the normal (radial) field 
components, as already noted before. 

The above expressions for a and c have appar- 
ent singularities at the poles, where they actually 
reduce in the generic case to resolved indetermina- 
cies with ^{(j), 9) and 6((/), 9) proportional to sm.9. 
Such indeterminacies are due to pole singularities 
of the global spherical system of coordinates rather 
than special properties of the field line mapping. 
Computationally, however, this requires care in 
evaluating Q near the poles, for example, by ex- 
cluding them from numerical grid. A more accu- 
rate approach is to use a covariant definition of 
Q with two overlapping regular coordinate charts 
covering the spherical boundaries (Titov 2007). 

5. FIELD STRUCTURE 

The photospheric distributions of Q and K 
are a two-dimensional imprint of the three- 
dimensional (3D) coronal field structure and so 
are very helpful for its analysis. Since their ex- 
pressions (18)-(23) are applicable to arbitrary field 
configurations, we have computed representative 
examples of the Q and K distributions for each of 
the three stages of magnetic field described above. 
In general, if the magnetic configuration has a 
separatrix surface (SS) due to the presence of an 
NP in the corona or a BP at the PIL, the field line 
mapping experiences a jump at the SS footprints 
and the corresponding Jacobian matrix is not de- 
fined there. Numerical derivatives, however, are 
estimated as a ratio of finite coordinates differ- 
ences, so their absolute values arc limited. Thus, 
numerically, the Q factor is always defined at the 
SS footprints and its narrow spiky distribution 
there helps identify NP SSs or BP SSs, if present. 
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Fig. 5. — Distributions of squashing factor Q (a,b) and expansion-contraction factor K (c) at the HFT 
footprints. On top of the grey shaded Q distribution (b) and red-blue shaded K distribution (c), the 
polarity inversion line (PIL, green line), the contours of Srl^^flp — 56i G, i = —3,..., 4 (blue and red 
lines for negative and positive values, respectively), two points of local minimum of B\^^j^^ (orange crosses) 
and footpoints of quasi-separator (magenta crosses) are superimposed. Only the parts of the distributions 
corresponding to the closed field lines are shown. 
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On the other hand, large values of Q, smoothly 
distributed on the photosphere in narrow strips, 
correspond to the footprints of QSLs (Titov et al. 
2002). These general properties of the Q distribu- 
tion are fully confirmed in the structural analysis 
of our field model. 

5.1. The initial potential field 

As expected from Section 2, the structure of 
the initial potential field in the AR turns out to 
be rather similar for both our models. Therefore 
we will describe further only the structure whose 
background field is determined in the framework 
of the PFSS model. The computed Q distribu- 
tion reveals two very sharp ridges with maxima 

lO^'', one for each magnetic polarity (see Fig- 
ure 5a, b). Each ridge spreads from the respec- 
tive bipole spot toward the pole region, curving in 
a sickle-like shape. They trace the footprints of 
the QSLs which combine into one flux tube, called 
an HFT. The cross-section of the HFT varies in 
a rather nontrivial way: neglecting some curving 
and twisting of the HFT, its cross-section can be 
represented as follows 

(24, 

This is possible because the cross-section shrinks 
in the longitudinal direction and stretches in the 
transverse direction when following the fleld lines. 
The dependence of such a variation on the arc 
length of the field lines has an exponential charac- 
ter. 

Because of this, resolving the internal structure 
of our HFT is numerically challenging. For this 
purpose, a curvilinear and nonuniform adaptive 
mesh with a minimal cell-size of < lO~^i?0 has 
been developed. To control the effect of trunca- 
tion errors, the computation of Q has been per- 
formed both with 10 and 12 significant digits by 
using the computer algebraic system Maple. The 
resolved profile of Q across the ridge of the Q dis- 
tribution in its highest part is shown in Figure 6. 
Similar to the case of a quadrupole magnetic field 
(Aulanier et al. 2005), the Q distribution has two 
peaks, which are probably related to the presence 
of two very localized minima of |B| at the HFT 
footprints (see below). Therefore, strictly speak- 
ing, the Q distribution has two sharp ridges closely 
bordering one another in each of the strips. How- 
ever, because this fine-scale feature of the HFT 



striicture is not significant for the results we report 
here, we neglect it in the subsequent discussion. 
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Fig. 6. — Profile of the logQ distribution in 
the perpendicular direction to the HFT footprint; 
s is the distance in this direction from one of 
the footpoints of the quasi-separator (</>, 9) ~ 
(3.251,0.935) (the north-east magenta cross in 
Figures 5 (b,c)). 

There is also a strong similarity between K dis- 
tributions of the quadrupole (Titov et al. 2002) 
and our configuration. In both cases, a very steep 
gradient of K exists across the ridges of the Q dis- 
tribution. On both sides of the ridges in the mid- 
dle part of the strips, there are areas of negative K 
(blue colored in Figure 5c). If crossing the strips 
in these areas outward from the PIL, the K dis- 
tribution with moderate negative values steepen 
first into a sharp front of highly negative values. 
It is followed then by a sharper front of positive 
values (red colored but not resolved in Figure 5c), 
which finally drops back to moderate negative val- 
ues. If moving along the strips from their middle 
to their ends, the areas and the fronts of negative 
K gradually disappear, while the fronts of posi- 
tive K increase in amplitude and broaden up to 
the entire width of the strips. 

By definition, positive and negative values of K 
correspond to expansion and contraction, respec- 
tively, of elemental flux tubes. This behavior of 
the flux tubes is opposite at the conjugate foot- 
points: if a given flux tube expands toward one 
footpoint, it contracts toward the other. Thus, 
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Fig. 7. — Map of the superimposed grey and red-blue shaded, respectively, Q and K distributions (see 
Figures 5b and 5c) on the solar surface (a); the process of assembling of the quasi-separatrix surfaces (QSSs) 
from the field lines launched at two ridges (b,c) of the Q distribution to give the entire structure (d). The 
intersection of the QSSs determines a quasi-separator field line (magenta solid (c) or dashed (b,d) line) with 
the footpoints marked by magenta crosses. Orange crosses are local minima of the photospheric \B\, and 
the green line is the PIL. 
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the conjugate footpoints of the field hnes are col- 
ored in Figure 5c in red and blue of equal intensity. 
In the Q distribution (Figure 5b), the conjugate 
footpoints are colored in grey of equal intensity. 
In combination with the PIL and contours of the 
photospheric Br superimposed on such distribu- 
tions, this provides a comprehensive graphical rep- 
resentation of the field line connectivity. One can 
see from this representation that most of the field 
lines in the HFT connect the end and mid parts of 
the strips or, in other words, the areas of stronger 
and weaker magnetic field. The field lines are or- 
ganized in the HFT in such a way that changing 
footpoints across the middle part of a given ridge 
causes swiping of the conjugate footpoint along 
the other ridge between the ends of the strips. In 
the case of the quadrupole configuration in each 
polarity, the ridges connect two magnetic flux con- 
centrations. In our present model, they similarly 
connect the flux concentrations of the AR with the 
respective polar magnetic "caps" . 

The grey-scaled Q distribution and the red- 
blue scaled K distribution can be superimposed 
on ((/), 6')-plane by using logical multiplication of 
pixel colors. Both distributions can then be com- 
bined in one colored graph, allowing one to easily 
compare these different characteristics of the field 
line mapping. To more clearly illustrate the prop- 
erties of Q and K, we have mapped the resulting 
image on the spherical surface and superimposed 
several contours of the corresponding photospheric 
Br component on it. The resulting image shown 
in Figure 7 summarizes the above discussed prop- 
erties of the field line mapping in the HFT. 

To better understand the HFT structure, it is 
also very instructive to trace the field lines start- 
ing at the ridges of the Q distribution. In ac- 
cordance with the above discussion, the conjugate 
footpoints are located not along the other ridge 
but rather across it and concentrated on a short 
arc passing through the point of a local maximum 
in Q. Each of such two sets of field lines forms 
a mid-surface of the corresponding QSL, which 
is further called a quasi-separatrix surface (QSS). 
These two surfaces intersect along a field line that 
we label the quasi-separator (Figure 7) for the fol- 
lowing two reasons. First, this field line is char- 
acterized by a maximum value of the squashing 
factor Q and it becomes a genuine separator field 
line in quadrupole configurations when their flux 



concentrations shrink to point sources (Titov et al. 
2002). Second, it is characteristic of the separator 
to pass through two magnetic null points of the 
configuration. The quasi-separator partly inherits 
this property in the sense that its footpoints are 
located near very localized minima of the magnetic 
field (Figure 5). 

In general, the presence of such minima in a 
given configuration is a good indication that the 
configuration also contains an HFT. Determining 
the minima of the photospheric distribution of |B| 
is much easier than analyzing the magnetic con- 
nectivity in the above manner. Therefore, an ini- 
tial scanning of the photospheric field minima can 
provide a quick estimate of the presence of HFTs 
in more complex configurations. It is necessary, 
however, to check additionally that the minima 
found are not due to the coronal magnetic nulls 
located right above the minima. This test can be 
done, for example, by plotting several iso-surfaces 
of \B\ = const > i?min to prove whether they do 
cover the minima and the gradient of |B| has the 
proper direction. In the above considered config- 
uration, such a test indeed confirms the absence 
of the coronal nulls above the photospheric field 
minima. 

The difference between the quadrupole case 
(Titov et al. 2002) and our model is mainly in the 
type of magnetic fiux distribution in one of the two 
bipoles of the configurations. In the quadrupole 
configuration, the flux is concentrated in two spots 
of the same AR, while in our model, it is dis- 
tributed over large solar hemispheres. Such a large 
length scale leads to a significant increase in our 
model of maximum values of Q, which are four or- 
ders of magnitude higher than in the quadrupole 
field. In all other respects, the magnetic fluxes in 
both conflgurations are partitioned by HFTs in a 
rather similar way. 

5.2. Sheared/twisted magnetic field 

Figure 8 (a) (c) shows the photospheric Q dis- 
tribution together with some of the new QSSs 
developed in the configuration after it is subject 
to the shearing/ twisting deformation described in 
Section 3.1. It can be clearly seen from Figure 8 
that the Q distribution outside the AR changes 
very little after shearing, which implies a corre- 
spondingly small change in the large-scale HFT 
initially present in the potential field. However, 
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(c) 











Fig. 8. — Distributions of Q at the end of the shear/twist (a) and flux-canceUation (d) phases; the QSSs in 
the AR in the first (b) and second (e) cases, respectively, and the corresponding views (c,f) of the Q and Br 
distributions in the vicinity of the AR; the PIL is a thick green hne, the contours of Br\^^j^^ = 44z G with 
i = —4, . . . , 4 (a-c) and i = —2, . . . , 3 (d-f ) are shown in hght blue and pink colors for negative and positive 
values, respectively; the thick white line in (e,f) is the footprint of the BP SSs wrapping around the flux 
rope depicted on (e) by a fraction of field hues. 
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on a smaller length-scale inside the AR, the mag- 
netic structure changes radically: two new HFTs 
form at the border of the region that is most 
affected by the photosphcric twisting motions. 
These HFTs have modest maxima of Q (~ 10^) 
compared to the initial large-scale HFT, but they 
are much shorter and located in the region of 
largest magnetic field strength. Therefore these 
newly formed small-scale HFTs play a more im- 
portant role in the flux-canccllation phase than 
the initial large-scale HFT. Such an expectation is 
fully confirmed by the analysis of magnetic struc- 
ture in the flux-cancellation phase (see sections 5.3 
and 6). 

The formation of HFTs by a twisting motion is a 
purely 3D effect that has no analogue in the two- 
and-half dimensional case of translationally invari- 
ant flux tubes. 

The vortex motions defined by equation (10) and 

applied to opposite flux spots of the AR have dif- 
ferent properties at the conjugate footprints of the 
twisting flux tube. In particular, these vortices are 
not concentric, in the sense that their centers of 
rotation are not connected by a field line. In com- 
bination with other differences in the vortices, this 
causes not only twisting but also braiding of the 
field lines in the coronal volume above the fiux 
spots. Such a braiding interlocks a part of the 
field lines and forms the above mentioned small- 
scale HFTs with less than two windings of the field 
lines. 

In order to describe qualitatively this effect, we 
construct now a simple analytical example of the 
field line mapping produced by similar vortices. 
The initial field line mapping from positive to neg- 
ative polarity in the potential configuration can 
approximately be represented in a Cartesian pho- 
tospheric plane with the PIL at the 2/-axis by 

n.^(^)^(-;)^ P5) 

A circular vortex motion centered at some point 
''c = {xc^Vc) can be defined as a differential rota- 
tion 

^ f cos fit — sinfii \ / x — Xc \f2Q) 
y sin Vlt cos nt J y y — Uc J 
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Fig. 9. — Distributions of Q produced by concen- 
tric (a) and eccentric (b) vortices in the opposite 
photosphcric polarities at t = 4.0 with = ^ , 
where p is the distance from the centers ri and 
of the vortices. In both cases n = (1.05, 0), while 
r2 = (-1.05,0) in case (a), and r2 = (-1.05,0.4) 
in case (b). 
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acting on a flux spot for time t with angular ve- 
locity, say, 

n = exp {-4 [{x -Xc)^ + {y- y,f] ' } . (27) 

The size of the vortex 1.2) corresponds here 
to the size of the flux spots. For simplicity, let 
the vortices be the same in positive and negative 
polarities and centered at points ri and r2, re- 
spectively. Assuming now that the configuration 
evolves according to ideal MHD, so that the field 
line connectivity is preserved, we can compose our 
field line mapping at time t as 

Ht = 7^r2,f °^oo7^rl,-t. (28) 

This equation means that the conjugate footpoint 
is obtained in three steps: first, a given footpoint 
is tracked back in time by T^n, -t to get its initial 
location, from which Ho provides then the initial 
location of the conjugate footpoint, whose current 
location is finally produced by TZr2.t- 

Each of the three mappings TZr^^t, IIo, and 
TZri,-t are incompressible, and so is their com- 
position Ilf. This implies that |A| = 1 in equa- 
tion (18) and so the squashing factor Q coincides 
with A''^ defined by equation (19), where the el- 
ements of the Jacobian matrix can easily be ob- 
tained from the above definition of XIj. Figure 9 
represents two calculated examples of Q distribu- 
tions for symmetric and asymmetric cases, where 
r2 and no(ri) are equal and not equal, respec- 
tively. It is natural to call such cases correspond- 
ingly concentric and eccentric. One can see from 
these examples that even a relatively small "eccen- 
tricity" in the differential rotations at the opposite 
polarities causes a substantial effect. Compared 
to the concentric case, Q increases here by an or- 
der of magnitude when each core of the vortices 
makes approximately one-half turn. The larger 
is the eccentricity in the vortices, the stronger is 
the squashing of elemental flux tubes at certain 
layers of the resulting twisted configuration. This 
has important implications for the current density 
distribution in the configuration under study, and 
also in a more general context, which is further 
discussed in section 6. 

5.3. Structure in the flux-cancellation 
phase 

Figure 8 (d)-(f) shows the Q distribution and 
QSSs at the final phase of the computation, af- 



ter flux cancellation has occurred at the PIL. The 
small-scale HFTs that are formed after the shear- 
ing phase have now increased in size. This occurs 
partly because of a noticeable migration of the 
conjugate footprints outward from the PIL, which 
is accompanied by reconnection in the HFTs due 
to a large enhancement of the current density in 
layer-like structures that will be discussed in sec- 
tion 6. 

Our simulations confirm also that the pho- 
tospheric flux-canceling process builds a twisted 
magnetic flux rope at the PIL (van Ballegooijen 
& Martens 1989). The inverse- S-like footprint of 
the respective BP SSs appears at this phase in the 
photospheric Q distribution by outlining the foot- 
prints of the flux rope [Figure 8(e-f)] in the way 
predicted earlier for a similar twisted structure by 
Titov & Demoulin (1999). Here it is evident how 
universal our technique is it enables us to iden- 
tify and locate all the structural elements of the 
magnetic field under study, including HFTs and 
BP SSs, as well as flux ropes. As a matter of fact, 
the very definition of a flux rope in a given asym- 
metric 3D configuration can be made precise via 
this technique. The flux rope is no more than a 
fraction of twisted field lines delimited from the 
surrounding configuration by BP SSs and/or he- 
lical QSLs [see another example of the rope in 
(Titov 2007)]. The presence of such "dehmiters" 
in the field configuration makes it possible to con- 
sider this fraction of field lines as a separate object, 
namely, a flux rope. 

Unfortunately, the boundary conditions chosen 
for modeling the flux-canceling phase of this run 
did not produce an eruption, even after reducing 
the magnetic flux of the AR to half of its initial 
value. The fraction of the canceled flux that has 
been converted into the flux of the rope is not large 
enough to reach an eruptive unstable state in this 
particular computational run. So the evolution of 
the magnetic structure at the eruptive phase is still 
to be investigated, which we plan to study after 
clarifying the eruptive boundary conditions for the 
flux-canceling process. We expect that the evolu- 
tion of the field structure prior to the eruption will 
remain similar because the HFTs and BP SS we 
have found are stable structural features that will 
survive under a wide variety of conditions. 

One can anticipate also how the magnetic field 
structure will change during the initial stage of the 
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eruption process. Based on the results of Titov 
& Demoulin (1999), we expect that detachment 
of the unstable flux rope will cause first a bifur- 
cation of the BP SS into a pair of BP SSs in- 
tersecting along a generalized separator field line. 
These bifurcated BP SSs will subsequently trans- 
form into an HFT located below the rising flux 
rope, as demonstrated by Titov (2007) for a simi- 
lar equilibrium configuration. Such a transforma- 
tion, however, will be a highly dynamical process, 
pinching this new HFT into a vertical current layer 
(CL), similar to the one described in previous work 
(Amari ct al. 2000; KHcm ct al. 2004). The mag- 
netic reconnection in this CL ought to play a cru- 
cial role in making the event eruptive (see sections 
6 and 7). 

Another important conclusion that follows from 
our example is that potential and non-potential 
magnetic fields with the same pliotosphcric flux 
distributions may have very different structures in 
both a topological and geometrical sense. This 
example undermines therefore a simplifying as- 
sumption often used in the analysis of solar mag- 
netic fields, namely that the topologies calculated 
in potential and non-potential approximations are 
similar. Hudson & Wheatland (1999) also argued 
against this assumption by proposing a counterex- 
ample of a constant-a force-free field that has a 
qualitatively different topology compared to the 
potential field with the same photospheric nor- 
mal component of magnetic field. However, Brown 
& Priest (2000) rightly argued against this coun- 
terexample by pointing out that it corresponds 
to a symmetric configuration, which is topologi- 
cally unstable to small perturbations of the field. 
Their consideration of asymmetric example sup- 
ported the opposite conclusion that the potential 
field model generally provides a good approxima- 
tion for the topology of the observed configura- 
tions. Our configuration is also asymmetric, but 
it does not support this conclusion, at least for 
the AR, where it has essentially nonliner force- free 
field. This means that the results of topological 
analysis of realistic magnetic configurations based 
on the potential approximation can be misleading 
and should be handled with caution. 



6. FIELD LINE CONNECTIVITY AND 
CURRENT LAYERS 

Photospheric shearing and twisting motions in- 
fiate the magnetic field of the AR, which presses up 
against the overlying large-scale HFT. The HFT 
footprints, however, being mostly outside of the 
applied photospheric vortices, are essentially not 
affected by these motions, so the line-tied HFT 
is stretched and pinched to counter the magnetic 
pressure from the infiating AR field. This pro- 
cess implies an accumulation of the current den- 
sity in the HFT such that its X-type mid-cross- 
section is noticeably pinched in the vertical direc- 
tion into a )-(-like shape with a relatively strong 
CL at the center (Figure 10). The numerical grid 
used in our MHD code does not enable us to re- 
solve the resulting thin CL and its four adjoining 
branches in its neighborhood. Although the cell 
size of the grid in this region 5 x 10~^ Rq is sev- 
eral times smaller than the apparent thickness of 
the CL, it is still much larger than the charac- 
teristic thickness of the corresponding QSL (see 
Figure 6), estimated as Rq / VQmax ~ 10~^7?q 
(Titov et al. 2002). Nevertheless, the asymptotic 
tail distribution of such a CL is consistently re- 
trieved by the code with an evident correlation 
between local maxima of squashing factor Q and 
parameter a = j\\/B at both shearing/twisting 
and flux-canceling phases. A similar correlation 
has been found by Aulanier et al. (2005) in MHD 
simulations of the quadrupole magnetic configu- 
ration in Cartesian geometry, where the corre- 
sponding HFT has been perturbed by distorting 
its footprints rather than inflating the underly- 
ing structure as in our case. This is consistent 
with the structural similarity of our model and 
the quadrupole configuration as discussed in sec- 
tion 5.1. 

One can expect therefore that this CL would 
lead to reconnection if eruptive conditions were 
met at the cancellation phase, to produce a rapidly 
rising flux rope that would push up more strongly 
against the overlying field structure. Such a 
CL is an essential element of eruption models in 
quadrupole configurations with an emerging inner 
bipole without shear (Syrovatskii 1982) or with 
strong shear (Antiochos et al. 1999), or with an 
embedded flux rope (Mackay & van Ballegooijen 
2006a, b), as in our case. We will argue in section 
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Fig. 10. — The distribution of a = j\\/B (normalized to 1/Rq) in the middle cross-sectional plane of 
the large-scale QSSs (depicted by several field lines) demonstrates an accumulation of the current density 
perturbations in the large-scale HFT at the end of shear/twist phase. The corresponding photospheric 
distribution of Q and the PIL are the same as in Figure 8a. 
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7 that the rcconnection in this CL imposes an im- 
portant constraint on the eruption process in the 
event under study. 

As discussed in section 5.2, two small-scale 
HFTs are formed in the AR during the shear- 
ing/twisting phase and they survive during the 
flux-canccUation phase of our simulations [pan- 
els (a)-(c) and (d)-(f) of Figure 8 for the former 
and latter phases, respectively]. In both phases, 
these HFTs carry strong layer-like currents, as ev- 
idenced by the sharp angles at which the corre- 
sponding QSSs intersect each other (sec Figures 8b 
and Be). In current-free HFTs such angles would 
be close to right angles. Figure 11 shows an ev- 
ident correlation between local maxima of Q and 
a and therefore confirms that thin CLs are indeed 
formed along the HFTs. One can also see from 
this figure that a is positive in the CLs and that 
they are located at the border of the twisted AR 
field, with negative a in the bulk volume. So these 
CLs carry a significant part of the return current 
that shields the twisted AR field from a nearly po- 
tential background field. Their formation occurs 
in the coronal volume due to the interlocking of 
the field lines by a twisting pair of shearing mo- 
tions eccentrically applied to the conjugate foot- 
points. In essence, this is another natural real- 
ization of Parker's mechanism of current layer for- 
mations discussed by Titov et al. (2003) and Gals- 
gaard et al. (2003) for quadrupole configurations. 
This process has been called magnetic pinching. 

Comparing the shearing/twisting and flux- 
cancellation phases [panels (a)-(f) and (g)-(n), 
respectively, in Figure 11], we sec that the CLs 
increase in size in the latter phase, whcrc^ a no- 
ticeable migration of their footprints outward from 
the PIL takes place. Since the magnetic Reynolds 
number based on the local velocities and length 
scale is not smaller than ~ 50, this migration is 
likely because of the reconnection that occurs in 
the CLs rather than due to magnetic diffusion at 
the photospheric boundary. Such CLs help re- 
distribute magnetic field and current within the 
AR at the flux-cancellation phase. More detailed 
investigation of this process is required, especially 
for the case when eruption occurs. 

It is clear, however, that its significance ex- 
tends far beyond this particular simulation. In- 
deed, magnetic sunspots with twisting-type mo- 
tions are quite natural and are often observed on 



the Sun (Brown et al. 2003) . For this reason, such 
motions have been considered for a long time as 
an in situ mechanism to build up magnetic en- 
ergy in the solar corona. Our example suggests, 
however, that the differential rotation of sunspots 
may also release this energy in the accompanying 
CLs by magnetic rcconnection triggered when the 
pinching of HFTs reaches a critical level. 

It should also be noted that the largest values 
of a are obtained in our simulations along the flux 
rope throughout its cross-section (panels (g)-(n) 
in Figure 11) rather than in the CLs. There are 
two reasons for this: first, the perturbation of the 
initial configuration by shearing/ twisting motions 
may not be sufficiently high to form stronger CLs, 
and second, these CLs are not well-resolved in our 
simulations. The latter may also be the reason 
why a CL is not formed along the BP SS that 
wraps around the flux rope, although the forma- 
tion of such a CL is generally expected in this con- 
figuration (Titov & Demoulin 1999; Low & Berger 
2003), even if it evolves quasi-statically. This is 
consistent with the numerical simulations of Fan 
& Gibson (2004), where such a CL is formed along 
the BP SS only in the dynamical kink-unstable 
phase of the evolution of a similar configuration. 
An additional separate study is required to clarify 
these issues. 

7. SUMMARY AND DISCUSSION 

We have developed a simplified model of the 
coronal magnetic field prior to the CME on May 12 

1997 in the following three steps. First, this field 
is constructed by superimposing a global back- 
ground field Bq and a localized bipolar field Bar 
to model the active region (AR) in the potential 
approximation. The background field Bq is de- 
termined from the observed photospheric normal 
field averaged over the longitude of the Sun. The 
AR field is modeled by a subphotospheric dipole 
whose parameters were optimized to fit the mag- 
netic field obtained from an MDI magnetogram. 
To verify the impact of the solar wind on the field 
structure, Bq is calculated for two different sets 
of boundary conditions. In the first case, Bq is 
modeled as a closed potential field vanishing at in- 
finity. In the second case, Bq is constrained to be 
strictly radial and therefore open at the sphere of 
radius r = 2.5 Rq. The resulting field structures 
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Fig. 11. — Distributions of logQ [(a), (d), (g), (1)] vs. a = j\\/B [(b), (e), (h), (m)] in the vicinity of flux 
concentrations at the end of shear/twist [(a)-(f)] and flux-cancellation [(g)-(n)] phases (a is normahzed to 
1/Rq); panels (a)~-(c) and (g)-(k) show these distributions at the photospheric ((/), 0)-plane, while panels 
(d)-(f) and (l)-(k) show them at the (6',r)-plane aX 4> 3.017 (yellow dashed line in the photospheric 
distributions). Superposition of the corresponding semi-transparent logQ distributions and a distributions 
is depicted by panels (c), (f), (k), (n). The PlLs and contours of -Brlr^flg the photospheric distributions 
are the same as in Figure 8. 
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in the AR and its neighborhood arc quite simi- 
lar for both choices of boundary conditions, there- 
fore only the first (closed) potential configuration 
has been used in the next step of MHD modeling. 
In this second step, the potential field is quasi- 
statically sheared by photospheric vortex motions 
applied to the two flux concentrations of the AR. 
Finally, the resulting force-free field is evolved by 
canceling the photospheric flux with the help of 
an appropriate tangential electric field applied to 
the central part of the AR. 

To investigate the structure of the modeled con- 
figurations, the field line mapping technique has 
been generalized to the case of spherical geometry. 
In our potential field configuration, this technique 
reveals a largc-scak; hyperbolic flux tube (HFT) 
characterized by very large values of the squash- 
ing factor. As in quadrupole-type configurations, 
such an HFT consists of two intersecting quasi- 
separatrix layers (QSLs) that provide an approx- 
imate partition of magnetic fluxes between four 
photospheric field sources. In our case, the role of 
such sources is played by two flux concentrations 
of the AR and two photospheric polar regions. 

It is important that no other structural fea- 
tures, such as magnetic null points or "bald 
patches" , are found in the AR and its surrounding 
area of the initial potential configuration. In this 
respect, our initial configuration drastically differs 
from that used in the breakout model (Antiochos 
et al. 1999), where the null points and separator 
are key topological elements of the structure. On 
the other hand, our large-scale HFT is a gener- 
alized analogue of the separator field line (Titov 
et al. 2002) and is therefore^ a favorable site for 
current layer (CL) formation (Titov et al. 2003; 
Galsgaard et al. 2003) and magnetic rcconncc- 
tion (Hornig & Priest 2003: Aulanier et al. 2005; 
Pontin et al. 2005). The accmmilation of current 
density observed in the vicinity of this HFT dur- 
ing the shearing/ twisting phase of our simulations 
is in agreement with this point of view. The mid 
X-type cross-section of the HFT pinches in the 
process of this accumulation into a )-{-like shape 
with a strong CL at the center. 

This structure is preserved during the cancella- 
tion phase when the flux rope is formed. Had the 
eruption conditions been met in this phase, the 
flux rope that would form would most likely de- 
tach from the dense photosphere and rise by push- 



ing the overlying flcld upward. This would cause 
an additional pinching of the CL to subsequently 
trigger reconnection within it. Such a reconnect- 
ing CL is a basic clement of the breakout model 
(Antiochos et al. 1999) and we have demonstrated 
that it appears in the flux-cancellation model as 
well if the interaction of the AR with the global 
background field is included. With this general- 
ization, the only difference between the models is 
in their interpretation of the drivers of the erup- 
tion. The breakout model assumes that the driver 
is reconnection in the CL, allowing for an explo- 
sive expansion of the sheared arcade of the AR. 
The flux-cancellation model assumes instead that 
the driver is a rising flux rope unleashed from the 
strapping field of the AR. While the reconnecting 
CL is important in our generalized model, the lat- 
ter mechanism appears to be more likely to drive 
the eruption. 

The following estimate of the magnetic fluxes 
in the interacting parts of our configuration sug- 
gests a picture of how this occurs. The flux as- 
sociated with the background field is eight times 
larger than the initial flux in our AR. However, 
only certain fractions of these fluxes interact with 
each other. These fractions can be estimated from 
the geometrical properties of the quasi-separator 
shown in Figure 7. Since the longitudinal exten- 
sion of the quasi-separator is about 0.55 radians 
(Fig. 5b), the fraction of the background mag- 
netic flux overlying the quasi-separator is equal to 
8 X 0.55/(27r) « 0.7. We must take into account 
also that only field components perpendicular to 
the quasi-separator can participate in reconnec- 
tion (reducing the fraction of the overlying flux 
from 0.7 to 0.4), and that the flux of the AR is 
reduced to half of its initial value at the end of 
the flux-cancellation phase. With these correc- 
tions, we obtain the interacting fractions of the 
background and AR fluxes expressed by the ratio 
$i0 : $iAR = 0.4 : 0.5. In other words, the in- 
teracting fluxes of the background and AR field 
are comparable. It implies that a substantial part 
of $iAR must be converted into the magnetic flux 
associated with the rising unstable flux rope in 
order to allow it to break through the overlying 
background field. Such a conversion is likely pro- 
vided by magnetic reconnection in the "vertical" 
CL below the rope, as observed in previous sim- 
ulations by Amari et al. (2000) and Kliem et al. 
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(2004). This process ought to replenish the mag- 
netic flux in the rope that is gradually destroyed 
by reconnection in the "horizontal" CL above the 
rising rope. 

We have also demonstrated that two small-scale 
HFTs pinched into thin CLs are formed by nonuni- 
form vortex motions applied to the flux spots of 
the AR. The reconnection in these CLs helps re- 
distribute the magnetic flux between the AR and 
the surrounding field during the flux-cancellation 
phase. Other possible implications of their pres- 
ence should be explored in a more accurate model 
of the magnetic field evolution in this event. How- 
ever, beyond its manifestation in this particular 
event, the discovered process of current layer for- 
mation by rotation of the flux spots has a broader 
significance. This process seems to be very im- 
portant for understanding other eruptive/flaring 
events driven by rotating sunspots (Brown et al. 
2003). Our example shows that their rotation can 
serve as a mechanism of both the storage of free 
magnetic energy and its release via gradual forma- 
tion of the above CLs in the solar corona. More 
detailed studies are certainly required in this di- 
rection. 

We have also confirmed that the photospheric 
flux-canceling process builds a magnetic flux rope 
at the polarity inversion line (van Ballegooijen & 
Martens 1989). This flux rope appears as a well- 
defined element of the structure owing to the bald 
patch separatrix surface that wraps around the 
rope. The magnetic structure near the rope is 
topologically equivalent to the twisted configura- 
tion described by Titov & Demoulin (1999). In 
particular, the bald patch separatrix surface has 
a characteristic inversc-S-like shape. This type of 
sigmoid corresponds to the negative sign of the 
current helicity in the flux rope, which is in accor- 
dance with observations (Liu 2004; Pevtsov et al. 
1997). 

V.S. Titov is thankful to Bernhard Kliem for 
several useful discussions on magnetic eruptions. 
This research was supported by NASA and the 
Center for Integrated Space Weather Modeling (an 
NSF Science and Technology Center). 
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A. PARAMETERS OF THE FIELD MODEL 



Tabic 1 prcsciiits the; cocifficicints of spherical harmonics for potential background magnetic fields Bq with 
two different types of boundary conditions. In both cases i?r© is prescribed at r = -R© = 1 by defining the 
Neumann boundary condition (4) for the harmonic potential Fq. In the first case, the additional asymptotic 
condition at infinity Bq\^^^ ^ is required to provide a standard exterior Neumann problem for the 
Laplace equation with the spherical boundary r = Rq. In the second case, the Dirichlet boundary condition 
^Q\r=R = at r = = 2.5 is imposed instead to define a mixed boundary value problem for Fq in the 
coronal volume Rq < r < i?ss, which is a so-called PFSS model. 

Table 2 presents the optimized parameters of the magnetic dipole in spherical coordinates to model the 
active-region field Bar through the gradient of the respective dipole potential (9). 
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Table 1 

Coefficients of spherical harmonics for two types of Fq 



n 






{R^, = 2.5Rq) 








ct 




1 


1.92227302 


1.86264984 


-0.119209589 


2 


-2.690655399 x 10"^ 


-2.672416196 x 10"^ 


2.73655421 x 10 


-4 


3 


1.33104612 


1.32943307 


-2.17814313 x 10" 


-3 


4 


-4.470680688 x 10"^ 


-4.46974778 x 10"^ 


1.171717562 X 10 


-5 


5 


0.480835203 


0.480804328 


-2.01663951 X 10 


-5 


6 


-2.583133904 X 10"^ 


-2.583117439 x 10~2 


1.733500769 x 10 


-7 


7 


1.1998345 X 10-2 


1.1982269 X 10-2 


-1.2865865 x 10 


-8 


8 


-5.299226134 x 10"^ 


-5.299194276 x 10"^ 


9.103946449 x 10" 


10 


9 


-7.47279826 x 10~2 


-7.47140767 x 10"^ 


2.0537249 x 10" 


-9 


10 


-1.633068626 x 10"^ 


-1.540580563 x 10"^ 


6.775544972 x 10" 


16 


11 


-2.40520137 x 10"^ 


-2.404263019 x 10"^ 


1.691849693 x 10" 


11 



Table 2 

Optimized parameters of the dipole 













0.9598074306 




0.9595343404 


0d 


1.186114672 




1.185997083 


4>d 


3.025771372 




3.025680793 


rrir 


1.740210202 X 10' 


-3 


1.705788497 x lO'^ 


me 


5.831919762 x 10 


-3 


5.84995686 x lO'^ 




1.800868961 x 10" 


-2 


1.811699206 X 10"^ 
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